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Abstract 

We present an undated comprehensive analysis for the simplest dark matter model in which a real singlet 
scalar with a Z2 symmetry is introduced to extend the standard model. According to the observed dark 
matter abundance, we predict the dark matter direct and indirect detection cross sections for the whole 
parameter space. The Breit-Wigner resonance effect has been considered to calculate the thermally averaged 
annihilation cross section. It is found that three regions can be excluded by the current direct and indirect 
dark matter search experiments. In addition, we also discuss the implication of this model for the Higgs 
searches at colliders. 
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I. INTRODUCTION 



The existence of dark matter (DM) is by now well confirmed UJ]. The recent cosmological 
observations have helped to establish the concordance cosmological model where the present Uni- 
verse consists of about 73% dark energy, 23% dark matter, and 4% atoms M\- Currently, several 
DM search experiments^ have found possible DM signals. The indirect DM detection experiments 
PAMELA yO, Fermi [143 and ATIC [|5D have observed the cosmic electron anomalies which can 
be explained by annihilating or decaying DM models The direct DM detection experiment 
CDMS II bl] observed two possible events in the signal region. In addition, the DAMA/LIBRA 
[|8|] and CoGeNT [|91] results favor a light DM candidate with a very large DM-nuclei elastic scat- 
tering cross section. 

In the standard model (SM) of particle physics, there is no candidate for dark matter. Therefore, 
one has to extend the SM to account for the dark matter. It is well known that the simplest DM 
model can be constructed by adding a real singlet scalar S to the standard model. In this model, 
a discrete Z2 symmetry (5 -S) has to be introduced to make the DM candidate S stable. Al- 
though this model is very simple, it is phenomenologically interesting. Motivated by the simplicity 
and predictability, a number of authors have explored its phenomenology [10-1^. These research 
works are very helpful for us to understand some more complicated DM models. 

In this paper, we try to give a comprehensive analysis for the real singlet scalar DM model. In- 
stead of assuming the fixed Higgs mass, we scan the whole parameter space of the DM and Higgs 
masses. For the resonance case, one should consider the Breit-Wigner resonance effect. Once 
the coupling between DM particle and Higgs boson is derived from the observed DM abundance, 
the DM elastic scattering cross section on a nucleon and the DM annihilation cross section in the 
galactic halo can be straightly calculated. The current DM direct and indirect detection experi- 
ments can be used to constrain the model parameter space. In addition, we also discuss the Higgs 
searches at colliders when the Higgs can decay into two DM particles. This paper is organized as 
follows: In Sec. II, we outline the main features of this model and give the DM annihilation cross 
sections. In Sec. Ill, we investigate the DM direct search, the DM indirect search and the collider 
implications. Some discussions and conclusions are given in Sec. IV. 
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II. THE REAL SINGLET SCALAR DARK MATTER MODEL 



In the real singlet scalar DM model, the Lagrangian reads 



£ = £sM + -duS&'S - - - AS'WH , 
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where H is the SM Higgs doublet. The linear and cubic terms of the scalar S are forbidden by the 
Z2 symmetry 5^-5. Then S has a vanishing vacuum expectation value (VEV) which ensures 
the DM candidate S stable. As describes the DM self-interaction strength which is independent 
of the DM annihilation. The observations of galactic DM halos and the dynamics of the Bullet 
cluster may constrain As when the DM mass is the order of 1 to 100 MeV lll7n . Notice that the 
one-loop vacuum stability can lead to a lower bound on the DM mass for a given As nl5\\ . It is 
clear that the DM-Higgs coupling A is the only one free parameter to control the DM annihilation. 
After the spontaneous symmetry breaking, one can obtain the DM mass = + A v^^ with 
vew = 246 GeV. It is natural for us to take mo in the range of a few GeV and a few hundred GeV. 
In this paper, we focus on 10 GeV < mo < 200 GeV. In addition, the Higgs mass m/j is also an 
important parameter for our numerical calculation. The 95% confidence-level (CL) lower bound 
on the Higgs mass is m/, > 114.4 GeV given by the LEP accelerator lllSn . The upper limit is 
mh< 186 GeV when we consider the precision electroweak data and the LEP direct lower limit 
il9|] . Therefore, we choose 1 14.4 GeV < irih < 186 GeV. It is worthwhile to stress that the current 



Tevatron experiments CDF and DO have excluded 162 GeV < m/, < 166 GeV [|20h . 

The real singlet scalar DM model is very simple and has only three free parameters mo, m/, and 
A. Based on the DM mass m^, the DM candidate S may annihilate into fermion pairs, gauge boson 
pairs or Higgs pairs. The annihilation cross sections & = 4EiE2crv {Ei and E2 are the energies of 
two incoming DM particles) for different annihilation channels have the following forms: 
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where s is the squared center-of-mass energy. The quantity F is defined as F(^) = arctanh(^)/^ 
with ^ = ^{s - Am^{s - 4mj^)/(s-2ml). When the DM annihilation nears a resonance, we should 
know the Higgs decay width F/, which is given by 

I,mUml-4m]y-^ ml j Am^ ( 4ml \lm%\ 
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Here we have included the decay mode /i ^ 55 if m/, > Imo. 



ni. DARK MATTER SEARCHES 



A. Dark matter relic density 



The thermal- average of the annihilation cross section times the relative velocity (crv) is a key 
quantity in the determination of the DM cosmic relic abundance. We adopt the usual single-integral 



formula for (crv) [12111 : 



with 
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&(s) = ag^ (9) 

where Ki(x) and ^2(-'c) are the modified Bessel functions, x = mo IT and g, = 1 is the internal 
degrees of freedom for the scalar dark matter 5 . Using the annihilation cross section & in Eqs. 
dlHS]), one can numerically calculate the thermally averaged annihilation cross section (crv) by use 
of the above formulas. 



The evolution of the DM abundance is given by the following Boltzmann equation [|22h : 

dY X six) T T 

where Y = n/s(x) denotes the dark matter number density. The entropy density s{x) and the Hubble 
parameter H evaluated at a: = 1 are given by 

2_7r^m| 
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where A/pL - 1 .22 x 10'^ GeV is the Planck energy, is the total number of effectively relativistic 
degrees of freedom. The numerical results of g^, have been presented in Ref. [|23n . Here we take 
the QCD phase transition temperature to be 150 MeV. 
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FIG. 1: The predicted coupling A as a. function of the DM mass mo from the observed DM abundance for 
the Mfi = 120 GeV and nif, = 180 GeV cases. The dashed line denotes the constant (crv> case when its value 
is taken at the usual freeze-out temperature x = 20. 



Using the result Yq of the integration of Eq. (fTOl) . we may obtain the DM relic density Q.ph^ = 
2.74X 10^ Yq mo/GeV. In terms of the observed DM abundance 0.1088 < Qoh^ < 0.1 158 [2], one 
can calculate the DM-Higgs coupling A for the given and m/,. As shown in Fig. [H the observed 
DM abundance requires A ~ (9(10 - 10^'). It is well known that the annihilation cross section 
cr will become larger for the same coupling when the annihilation process nears a resonance. This 
feature implies that there is a very small coupling when 0.8 m/, < 2mj) < m/,. This region is 
named as the resonance region in the following parts of this paper. It should be mentioned that the 
thermally averaged annihilation cross section (crv) will significantly change as the evolution of the 
Universe when the DM particle is nearly one-half the mass of a resonance [|24|] . This is the Breit- 
Wigner resonance effect which has been used to explain the recent PAMELA, ATIC and Fermi 
anomalies. Here we have considered the Breit-Wigner resonance effect for the determination of 
the coupling A. If we neglect the Breit-Wigner resonance effect and only consider the resonance 



contribution, the predicted coupling A will has distinct differences with the previous results. For 
illustration, we also present the corresponding A in Fig. \T\hy use of the constant (crv) at x = 20. 
In this case, the predicted A may differ from the correct one by more than a factor of 3 for the 
resonance region. Since the DM direct and indirect detection cross sections are proportional to A^, 
one can derive the larger differences. 



B. Dark matter direct search 



For the scalar dark matter, the DM elastic scattering cross section on a nucleus N is spin- 
independent, which is given by yj] 

4M^(N) 

CTN = —(Zfp + (A- Z)Uf , (13) 

n 

where M(N) = moMjs{/{mj} + Mj^) and Mj^ is the target nucleus mass. Z and A-Z are the numbers 
of protons and neutrons in the nucleus, fp^, is the coupling between DM and protons or neutrons, 
given by 

fp,n = ^ frq 27-^'"'<^ ^ ^ ^'^~m~' ^^^^ 

q=u,d,s 1 q=c,b,t ^ 

where f^!^^ = 0.020 + 0.004, /^J = 0.026 + 0.005, /}^^ = 0.118 + 0.062, fi^ = 0.014 + 0.003, 



f^"J = 0.036 + 0.008 and = 0.1 18 + 0.062 UM- The coupling /^^'^ between DM and gluons 
from heavy quark loops is obtained from = 1 - 'Z.g^uM,^ fj'q"^^ which leads to fj^ ~ 0.84 and 
fj^ « 0.83. The results of DM-nucleus elastic scattering experiments are presented in the form 
of a normalized DM-nucleon scattering cross section erf ^ in the spin-independent case, which is 
straightforward 

/ \ 2 / \ 2 
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"^'^ \ niD + m„ ) A^M^N) * ;r \ + m„ j ^" ' ^^^^ 
where we have used fp « /„. 

In the real singlet scalar DM model, the DM candidate S interacts with nucleus N through 
^-channel Higgs exchange. In this case, the DM-quark coupling in Eq. (fT4l) is given by 

A m„ 

a, = . (16) 

niD ml 

Using the predicted A from the observed DM abundance, one can calculate the spin-independent 
DM-nucleon elastic scattering cross section erf ^ for the given mo and m/,. We perform a numerical 




FIG. 2: The predicted DM-nucleon scattering cross section erf/ for the whole parameter space. The short 
dotted lines with arrowhead indicate the excluded regions from the DM direct search experiments CDMS II 
and XENON 10. The region among two dashed lines (162 GeV < m/, < 166 GeV) can be excluded by the 
Tevatron experiments CDF and DO. Two purple short dashed lines describe the minimum mo allowed by 
both the DM observed abundance and the vacuum stability/pertubativity for As = 10"^ and As = 0.4. The 
right panel corresponds to two fixed Higgs mass cases with current and future experimental upper bounds. 

scan over the whole parameter space of and m/,. The numerical results are shown in Fig. [2l 
For illustration, we have also plotted erf/ as a function of mo for the m/j = 120 GeV and mh_= 180 
GeV cases. In view of the current DM direct detection upper bounds from CDMS II [Tl] and 
XENON 10 [|26n . we find that two regions indicated by short dotted lines can be excluded. The 



future experiments XENONIOO IZZD, CDMS 100 kg and XENONIT can cover most 
parts of the allowed parameter space. For the resonance region, the predicted erf/ is far below 
the current experimental upper bounds. In Fig. [2l we also plot the minimum mo allowed by 
both the DM observed abundance and the vacuum stability/pertubativity for the DM self-coupling 
As = 10"^ and As = 0.4 when the cutoff scale is taken to be 1 TeV [llSll . Increasing the cutoff 
scale, one can derive much stronger bounds lllSn . In addition, the region among two dashed lines 



(162 GeV < m/, < 166 GeV) can be excluded by the Tevatron experiments CDF and DO. 
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FIG. 3: The predicted thermally averaged DM annihilation cross section (crv>o- The right panel corresponds 
to two fixed Higgs mass cases. 

C. Dark matter indirect search 

The DM indirect search experiments are designed to detect the DM annihilation productions 
which include neutrinos, gamma rays, electrons, positrons, protons and antiprotons. Since these 
methods are complementary to the direct detection, it will be very important for us to investigate 
the DM indirect detection. In order to derive the correct relic density, we have calculated the ther- 
mally averaged annihilation cross section (crv) which depends on the temperature of the Universe. 
On the other hand, (crv) also determines the DM annihilation rate in the galactic halo. The only 
difference among the above two cases is the temperature T. For the relic density, (crv) is usually 
evaluated at the freeze-out temperature x « 20 (the averaged velocity v « yTiJx). The dark matter 
annihilation in the galactic halo occurs at x « 3 x 10^ (v « 10^^). Therefore we should calculate the 
thermally averaged annihilation cross section at x « 3 x 10^, namely (crv)o. The numerical results 
have been shown in Fig. [31 One may find 1 x 10"^^ cm^ sec"^ < (crv)o < 3 x 10^^^ cm^ sec"' for 
most parts of the parameter space, which is consistent with the usual 5- wave annihilation cross sec- 
tion (crv) « 3 X 10"^^ cm^ sec"' at the freeze-out temperature x ~ 20. However, the Breit-Wigner 
resonance effect can enhance or suppress (crv) for the resonance case [2^- As shown in Fig. 
|3l (crv) in the red and yellow regions is enhanced by the Breit-Wigner enhancement mechanism. 
Contrarily, (crv) is suppressed ((crv) <^ 1 x 10"^^ cm^ sec"') when double DM mass 2mj) is slightly 
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less than the Higgs mass m/,. In this case, it is very difficult for us to detect the signals of the DM 
annihilation. The vertical cyan region around 80 GeV has smaller (crv)o which originates from 
the threshold effect BOll . If is slightly less than the W boson mass, the channel SS ^ W'^W~ 
is open at high temperature. It dominates the total thermally averaged annihilation cross section 
and determines the DM relic density. However, this channel is forbidden in the galactic halo (low 
relative velocity). At this moment, the channel SS ^ bb has the dominant contribution. This 
feature can be well understood from Fig. [3] (right panel). 

In our model, the DM annihilation can generate primary antiprotons which can be detected 
by the DM indirect search experiments. Recently, the R^.MELA collaboration reports that the 
observed antiproton data is consistent with the usual estimation value of the secondary antiproton 
[|3n. Therefore one can use the R\MELA antiproton measurements to constrain (crv)o. In Fig. 
|3l we have also plotted the maximum allowed {<tv)o for the MIN, MED and MAX antiproton 
prop agation models [|l4ll . Notice that a fixed Higgs mass m/j = 120 GeV has been assumed in Ref. 
[|l4ll . The above upper bounds are still valid for our analysis when m^, < 1 14.4 GeV. Then we can 
find that a very narrow region can be excluded by the R\MELA antiproton data in our model. This 
feature is not shown in Ref. [|l4|]. In fact, the width of this excluded region is about 0.4 GeV if we 
require (crv)o < 10"^^ cm^ sec~^ For the MED and MAX propagation cases, the sensitivity of the 
future experiment AMS-02 [31] may reach <crv)o ~ <9(10"^^ - 10"^^) cm^ sec^' [Qj], which can 
cover most parts of the whole parameter space as shown in Fig. |3l 



D. Implications on the Higgs search 

Since the DM candidate S has substantial coupling to nucleons via Higgs boson exchange, 
they can be produced in high energy collider experiments such as the Tevatron and Large Hadron 
Collider (LHC). If /m/, > 2mo, the decay channel /z ^ 25 is kinematically allowed. In this case, 
two DM particle production may be associated with the Higgs production. The produced DM 
particles are invisible and have the missing energy signals. This will affect the usual SM Higgs 
searches at the Tevatron and LHC. To describe this effect, we calculate the branching ratio of the 



Higgs visible decay llllll 



BRvisible - — • (17) 
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FIG. 4: The predicted branching ratio of the Higgs visible decay BRyisibie- The right panel corresponds to 
two fixed Higgs mass cases. 

The numerical results have been shown in Fig. |4l Considering the constraints from the DM direct 
search experiments, we find BRvisiWe ^ 0-3. In fact, we have BRyisibie > 0.9 for most parts of 
the allowed parameter space. If the future CDMS 100 kg does not observe the DM signals, the 
region less than 0.9 can be excluded. Notice that the invisible Higgs decays in the allowed region 



0.3 < BRvisibie ^ 0.75 can be identified at the ATLAS with an integrated luminosity of 30 fb [|32|] 



The decay mode h ^ SS can lower the statistical significance of the traditional Higgs search at 



the CMS lllBn . A combined analysis of the traditional visible search and the invisible search at the 
LHC can enhance the discovery reach and the possibility of distinguishing this DM model from 
the SM 0. 

IV. DISCUSSION AND CONCLUSION 

We have detailedly discussed the 10 GeV < nio < 200 GeV case. As shown in Fig. [2l 
one can obtain the very large DM-nucleon elastic scattering cross section crfj ~ O(\0~^^cm^) for 
mo ~ 10 GeV. In this case, the light DM particle S can explain the DAM A/LIBRA and CoGeNT 
experiments [16]. It should be mentioned that this explanation is consistent with the XENON 10 



and CDMS null results ||331]. Currently, the XENON 100 preliminary results challenge the DM 



interpretation of the DAMA/LIBRA and CoGeNT signals 



However, there are bifurcations 
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on the choice of the ratio between electron equivalent energy and nuclear recoil energy ySQ. If 
niD > 200 GeV, one will not meet the resonance and the new annihilation channels. Then (crv)o « 
2.3 X 10"^^ cm^ sec"' and 10""*^ cm^ < erf/ < 10'"^ cm^ can be derived. 

In conclusion, we have made an undated comprehensive analysis for the whole parameter space 
of the real singlet scalar DM model. To satisfy the observed DM abundance, we numerically solve 
the Boltzmann equation and predict the DM-Higgs coupling A which determines the DM direct and 
indirect detection rates. We demonstrate that the Breit-Wigner resonance effect can significantly 
change the predicted coupling A by more than a factor of 3 for the resonance region. The spin- 
independent DM-nucleon elastic scattering cross section crf^' has been presented for the whole 
parameter space of nio and mf,. One may find that the current experimental upper bounds from 
CDMS II and XENON 10 can exclude two regions. For the DM indirect detection, we calculate 
the thermally averaged annihilation cross section (crv)o which can be enhanced or suppressed by 
the Breit-Wigner resonance effect. We find that a very narrow region can be excluded by the 
PAMELA antiproton measurements. One should notice that the predicted erf/ and (crv)o are very 
small for the resonance region. In this case, it is very difficult for us to detect the signals of the DM 
annihilation. However, we still have possibility to test the resonance region as detector masses of 
DM direct search experiments become larger. When the decay channel h ^ SS is kinematically 
allowed, we find that the allowed branching ratio of the Higgs visible decay satisfy BRvisibie <: 0.3. 
If the future CDMS 100 kg does not observe the DM signals, the region less than 0.9 can be 
excluded. 
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